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Abstract 



We connect the observed under-abundances of Li and Be in dwarfs, with recent 
results on nuclear cross sections at low energies: for collisions of protons with 
atomic or molecular targets, the measured cross sections seem too high with 
respect to extrapolations for bare nuclei. Phenomenologically, these anoma- 
lous nuclear interactions can be described in terms of an effective screening 
potential Ui a b in the range of few hundred eV: in the presence of the elec- 
tron cloud, nuclei become more transparent to each other as if the effective 
collision energy is aumented by £// a fc- This implies that fusion cross sections 
are enlarged and at the same time elastic cross sections are lowered. If some- 
thing similar occurs in stellar plasma, the nuclear burning temperatures are 
lowered, whereas diffusion processes are enhanced. We find that the observed 
Li and Be abundances in the Hyades and in the Sun can be reproduced for 
effective screening potentials of the plasma in the range of 600-700 eV, close 
to that found by experiments in the laboratory. 

subject headings: stars: abundance-nuclear reactions- diffusion. 
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1. INTRODUCTION 



This paper is an attempt to connect two anomalies, one in stars and the other in the 
laboratory: 

a) the abundances of Li and Be in dwarfs are generally too small with respect to the predic- 
tions of evolutionary calculations; 

b) the measured nuclear cross sections at low energies for collisions of protons with atomic 
(or molecular) targets, a at , are too large, with respect to extrapolations of data taken at 
higher energies, where one can assume to have collisions of protons with bare nuclei, cr nuc . 

Concerning a) the problem is to account for the observed depletion of light elements in 
photospheric regions, the temperature - in these regions and in the convective layers below - 
being too low for nuclear burning, according to the nuclear reaction rates adopted in current 
evolutionary codes. Two well known evidences of this situation are given by the Sun and 
the Hyades, see Michaud & Charbonneau (1991) for a review of observational data and 
theoretical approaches. 

We remark that the Li abundance problem is not alone, but has to be seen and explained 
in conjunction with the fate of other light elements, for example Be. The comparison between 
the present Sun and the meteoritic composition indicates that the solar photosphere is 
depleted by a factor about 100 for Li and about 2 for Be. 

Concerning b), the situation is that for several nuclear reactions, at collision energies 
E in the 10 — 100 KeV range the measured cross sections are too large with respect to 
expectation. If one writes 

Cat(E) = CT nuc (E + Ulab) , (1) 

one finds £7/ a & much larger than the values U e i given by theoretical calculations for the 
electron screening. For example, for collisions of protons with 6 Li and 7 Li nuclei one finds 
experimentally U[ a + b U 400 eV (Engstler et al 1992), whereas U^ +U 180 eV (Bracci et 
al. 1990). 

We call Uiab an "effective screening potential in the laboratory" and we consider it as a 
parameter accounting for the (presently) inadequate theoretical treatment of nuclear reac- 
tions at very low energies. We remark that the condition £7/ a & >> U e i is not an accident for 
just one reaction, but it seems to be a general problem, see e.g. Fiorentini, Kavanagh and 
Rolfs (1995). 

Experiments in the laboratory, if correct, indicate that some effects due to the electron 
cloud around the colliding nuclei has been neglected. May be that also in stars we are 
neglecting something, connected with the interactions of the nuclei with the plasma. The 
aim of this paper is to investigate such a possibility. 

As a working hypothesis, we thus make the following Ansatz: some anomaly also occurs 
in the stellar plasma, so that the nuclear cross sections in the plasma are given by 

Cpi(E) = cr nuc (E + U p i) (2) 

and we keep U p i, the "effective plasma screening potential" for the reaction, as a free variable 
allowing it to be in the range of a few hundred eV, i.e. significantly larger than found in 
standard treatments of plasma screening (as an example, for collisions of two nuclei with 
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charges Zi and Z 2 at the bottom of the solar convective region where the Debye length 
is Rr> 2. 10~ 8 cm, in the weak screening approximation one has U w = ZiZ 2 e 2 / Rr> 
7Z 1 Z 2 [eV]). 

We also allow U p i to be specific of the nuclear reaction, 

U f i = U*t Az , (3) 

but we assume that it is independent of plasma properties (at least in the region of interest 
to us), much in the same way as the observed £7/ a & values are approximately independent of 
the chemical or physical state of the target. 

Equation (2) implies that in stellar calculations the Maxwellian averaged burning rates 
at temperature T become: 

Ap,(T) = X nuc (T)exp(U pl /kT) (4) 

where X nuc are the burning rates for bare nuclei. This clearly gives significant changes in 
the temperature Tnb a t which nuclear burning becomes efficient (which we define from the 
condition A„„ c (T/v#)i0 = 1, t & being the Sun's age). For example, the burning temperature 
of 7 Li in the Sun, usually quoted as T/v#=2.5 (here and in the following T means the 
temperature in units of 10 6 K), becomes T/v#=2.1 for t/£ +L, =700 eV, whereas for 9 Be one 
moves from T/v#=3.5 to T/v#=3 passing from Upf Be =0 eV to Upf Be =700 eV. 

The basic question addressed in this paper is thus the following: can we find a range of 
U p i values which can account for the observational data in stars? 

Firstly, in section 2 we present a discussion of the observational data for 1M stars, in 
the conventional framework, with the following main conclusions. 

a) We exclude any attempt of explaining Be depletion in terms of mixing the photosphere 
with high temperature (T ss3.5) regions where Be can be burnt, since Li would be completely 
destroyed. 

b) As a consequence, Be has to be hidden (without being burnt) below the convective zone 
by diffusion processes. 

c) From the comparison essentially among meteorites, Hyades and the present Sun, we 
can estimate Li depletion in the pre-main sequence (PMS) and in the main sequence (MS) 
phases. Li depletion occurs in both phases. 

d) The same diffusive processes which hide Be are not sufficient to account for the Li deple- 
tion during the main sequence. Li has to be burnt also during the main sequence. 

In section 3 we summarize, from measurements in the laboratory, the knowledge of £7/ a & 
for the reactions of interest to us and we present some arguments for assuming that the 
plasma screening potential U p i is (roughly) independent of plasma properties, and similar 
(not necessarily equal) to £7/ a fc- 

We also show that introducing U p i significantly enhances the diffusion coefficients. The 
anomalously large measured nuclear cross sections seem to indicate that nuclei are more 
transparent to each other and this suggests that elastic cross sections are lowered, so that 
diffusion can be enhanced. 

In section 4 we investigate Li burning in the pre-main sequence and in the early main 
sequence phases. The study of a young open cluster such as the Hyades, for which many 
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observational data are available, is of particular interest since one can compare stars with dif- 
ferent masses while the other parameters (age, chemical composition...) are the same. Stars 
with different masses have different values of temperature at the bottom of the convective 
zone, stars with lower surface temperature being hotter at the bottom of the convective 
zone. In this way one can essentially explore the burning rate as a function of temperature. 
In a sense, this is the astrophysical equivalent of an experiment in the laboratory where 
the energy of the accelerated particles is varied so as to measure the energy dependence of 
nuclear cross sections. We will find that for Ufr Li « 600 - 700eV one can account for the 
observational data on the Hyades. 

In section 5 we discuss Li (and Be) depletion of the convective zone in the main sequence 
phase of the Sun. This depletion can be due to several processes: 

a) burning at the bottom of the convective zone, which clearly gets easier if the anomalous 
burning rates of equation (4) are considered; 

b) concentration driven diffusion: for Li, the presence of an anomalous screening potential 
results in a steep rise of the burning rate just below the convective zone; this induces a 
strong Li concentration gradient which can drive diffusion; 

c) gravitational and thermal diffusion. 

The efficiency of processes b) and c) depends on the values of the diffusion coefficients. 
With the help of the observed value for the Be depletion in the Sun, which we consistently 
assume to be due to diffusive processes, we will fix the diffusion coefficients for Li and we 
will study the effects of processes a)-c). We will show that again for U^ Ll ~600-700 eV we 
have a mechanism which can account for the observed solar Li depletion. 

Throughout all the discussion, the properties of the bottom of the convective region 
play a key role. Within a phenomenological approach, these properties will be determined 
from observational (helioseismoligical) data whenever possible, otherwise we will use results 
which are common to most solar model calculations. All this information is summarized 
in the Appendix B, whereas Appendix A contains a short outline of the evolutionary code 
FRANEC (Frascati Raphson Newton Evolutionary Code) we are using. 

In conclusion, the range of V^ L% which is needed for solar Li depletion in the main 
sequence overlaps with that reproducing the Li abundances in the Hyades and it is not far 
from the value of Uf^ Ll measured in the laboratory. A short discussion of future prospects 
is given in the final section. 

2. OBSERVATIONAL DATA ON LI AND BE DEPLETION IN 1M STARS 

Observational data on Li depletion in different astrophysical contexts have been recently 
reviewed by Michaud & Charbonneau (1991). From the available data we selected those 
pertaining to 1M stars in various stellar clusters. Assuming that the selected stars have 
similar chemical composition, the graph shown in Fig. 1 tells us the history of photospheric 
Li abundance in 1M star. The very stable structure of main sequence stars implies a 
constant efficiency of physical mechanisms and thus suggests a linear time dependence of 
photospheric abundances. Indeed, the data shown in Fig. 1, are well fitted by a straight 
line, 

yu{t) = (2.35 ± 0.2) - (0.27 ± 0.05)* , (5) 
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where as usual y a = 12 + log(N a /Nff), t is the age in Gy and the quoted errors are to be 
taken as indicative of the uncertainties of the observational data. There is a clear indication 
that some Litihum depletion occurs during the main sequence. For the solar age we get: 

Ay L *(MS) = 1.2 ± 0.25 (6) 

In Fig.l we have also shown the meteoritic value, which tentatively can be taken as repre- 
sentative of the initial composition. The linear fit, eq. (5) then implies that a significant Li 
depletion also occurred in the pre-main sequence phase: 

Ay L8 (PMS) = 1.0 ± 0.2 (7) 

We note that this result is consistent with observational data from TT tauri which for 1M© 
stars give only an upper bound Ay^ < 1 ) (Zappala 1972; Magazzu, Rebolo and Pavlenko 
1992). 

In stellar environments we do not have so far any observation about the relative abun- 
dances of 6 Li and 7 Li. The burning rates in the region of interest to us (T=2-4) can be 
espressed as (see Fig. 2): 

Kit 2 = Ba zP XT^z , (8) 

with: 

B* Li = 9.5 • 10~ r Gy' 1 ; ae Ll = 19.6 (9a) 



Br u = 1.2 • lO" 8 ^- 1 ; aru = 19.7, (9b) 

X being the Hydrogen mass fraction and p the density, here and in the following in g/cm 3 . 

The burning rate of 6 Li is thus about a factor hundred higher than that of 7 Li. The 
small initial 6 Li fraction (about 8% in meteorites), is thus completely destroyed when (and 
if) 7 Li is burnt. For this reason in the following we will neglect 6 Li and refer for brevity to 
Li when actually discussing 7 Li abundance. 

The few available data about Be in 1M stars are summarized in Table 1. Briefly, the 
solar abundance is one-halfof that in the meteorites, and data from Hyades for 1M stars 
cannot discriminate between the solar and the meteoritic values, at present. Anyhow, this is 
important to tell the fate of Be. In the temperature region of interest to us the conventional 
(U p i=0) burning rate can be expressed again using equation (8). When the two final channels 
a + 6 Li and d+2ct are summed one has: 

B, Be = 2.5 • lO" 13 ^- 1 ; a, Be = 23.5 (10) 

Thus the ratio between the burning rates of 7 Li and 9 Be is 

R = K + ul u /K + ul Be = 4-8 • io 4 t- 3 - 8 (11) 

Li and Be nuclei are close in mass and charge (also in any calculation the diffusion 
coefficients come out to be very similar for both nuclei) and it is natural to assume that 
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any mixing mechanism is the same for both nuclei, so that they will experience the same 
temperatures. It is then clear from equation (11), see also Fig. 2, that at temperatures such 
that Be can be burnt (T 3.5) Li would be completely destroyed!! 

This rules out any attempt to explain Be depletion in terms of mixing the photosphere 
with high temperature regions where Be can be burnt. 

There are only two possile outcomes: 

a) Be can be burnt at significantly lower temperatures, which are not too dangerous for Li 
survival. 

b) Be depletion is just a diffusion process: Be is hidden (not burnt) below the convective 
zone. 

As will be discussed in sections 4 and 5, Be could be burnt in 1M stars for U^ Be ~1700 
eV, a value which seems too high, although it is not excluded by the few experimental 
laboratory data presently available, see next section. We are thus left with the conclusion 
that Be is hidden just below the convective zone by diffusion. This requires the long times 
of the main sequence phase, and is due to pressure and thermal diffusion, concentration 
gradients being too low at the bottom of the convective zone in order to affect the Be 
diffusion. 

The same processes will also affect the fate of Li, due to the similarity of the two elements. 
We thus conclude that during the main sequence the Sun has lost a fraction 

Ay£ T (MS) « Ay Be = 0.27 ± 0.10 (12) 

where the index PT reminds the effects of pressure and temperature. The comparison 
between equation (6) and equation (12) shows that other effects are to be added, in order 
to account for the solar Li depletion during the MS phase. 

Before closing this section, let us summarize the phenomenological picture we presented, 
recalling also the results of current evolutionary calculations. 

i) There is observational evidence for a significant Li burning in the pre-main sequence phase 
of a 1M star. However, calculations in the standard framework (e.g. Proffitt & Michaud 
1989) do not achieve significant Li burning, essentially since the star is not hot enough 
during the short pre-main sequence times. Either the pre-main sequence models are wrong 
or Li can be burnt more easily than believed. 

ii) During the main sequence, Li diffusion due to PT effects, as estimated from data on 
Be, accounts just for a fraction of the Li depletion, and thus Li burning has to occur also 
in the main sequence. From helioseismology, and in agreement with several evolutionary 
calculations (see Appendix B), the temperature at the bottom of the convective zone of 
the Sun is determined to be in the range T{,=2. 1-2.3 and this temperature should not have 
changed by more than 10% during the previous Sun evolution. This leads again to the 
conclusion that L can be burnt at temperatures smaller than currently believed. 

iii) We ascribe Be depletion in the Sun to diffusion processes. As is well known (see section 
5 and Appendix B), theoretical calculations (neglecting turbulence) yield a too small Be 
depletion. 

A mechanism that at the same time would give lower burning temperatures and faster 
diffusion than in the standard framework would be welcome. 
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3. NUCLEAR BURNING IN THE LABORATORY: EXPERIMENTAL RESULTS 
AND SUGGESTIONS FOR STELLAR INTERIORS 



An extensive experimental investigation of nuclear reactions between charged particles 
has been performed at very low collision energies (EsslO-100 KeV) in the last few years, 
particularly by the Rolfs group (Assenbaum et al. 1987; Engstler et al. 1992 and references 
therein). 

The striking result is that at these low energies the nuclear cross sections are generally 
larger then expected by extrapolating data taken at higher energies, where the effect of the 
electron cloud surrounding the target nuclei can be neglected. By parametrizing the cross 
section as in equation (1), one extracts from the experiments the values reported in Table 2 
for the effective screening potential £7/ a fc- 

Electron screening effects have been calculated by using different approximations: adi- 
abatic/sudden (Bracci et al. 1990), classical trajectory Monte Carlo method (Bracci et 
al. 1989) and in a few cases with an ab initio quantum mechanical dynamical calculation 
(Bracci et al. 1991). The maximal value of the electron screening potential is obtained 
in the adiabatic limit, which corresponds to the maximal energy which, consistently with 
quantum mechanics, can be transferred from the electrons to the nuclear motion (Bracci et 
al 1989). In this limit, for collisions between two atomic systems A and B yielding at zero 
internuclear distance the compound atomic system A+B, the screening potential is given 
by: 

U^ +B = E(A + B)- E(A) - E(B) , (13) 

where E(i) are the atomic binding energies. 

Although the experimental errors are individually large, the measured values £7/ a & are 
systematically larger then U e \, see again Table 2. The reason for this discrepancy is not 
understood. We take the attitude that it is due to some inadequacy of the theoretical 
treatment and try to learn directly from experiments the properties of £7/ a fc- 

Reactions between H isotopes and Li isotopes are of particular interest for the present 
discussion. As one sees from Table 2, the values of U^ h +Ll are, whithin errors, i) independent 
of the isotopes which are interacting and ii) independent of the chemical state of the target. 
As remarked in Engstler et al. (1992), the first point suggests that the problem is not due to 
nuclear physics but it is related to interaction with the electron cloud. On the other hand, 
the second point shows that the effect is weakly sensitive to the detailed structure of the 
electron cloud. 

Concerning 9 Be(p,ai) 6 Li and 9 Be(p,d)2ai reactions, the only published data at low enough 
energies are from Sierk & Tombrello (1973). At the four lowest measured energies, data show 
an increase of the astrophysical factor corresponding to £7/ a & in the range of about 1 KeV, 
however the statistics is too poor for getting a definite answer. Preliminary results from 
Rolfs group (Zahnow et al. 1994) seem not to confirm the low energy enhancement. 

In a stellar plasma of interest to us (T=2-4, p (0.1 — 1), X ss0.7) only a small fraction of 
Li or Be atoms are not dissociated. On the other hand, it is worth observing that the plasma 
cloud around a Li nucleus is similar to the electron cloud around the target Li nucleus in the 
laboratory, as in both cases the mean interparticle distance is d 10~ 8 cm and the electron 
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velocity is v e i 10 8 cm/ sec. It is thus reasonable to assume an effective screening potential 
in the plasma V^ L% of the same order of magnitude as that measured in the laboratory, in 
the range - say - of a few hundred eV. 

In addition, d and v e i change weakly inside the plasma region defined above. Reminding 
the weak sensitivity of U^ Ll to the physical state of the target, it is thus natural to take 
Upi~ L * as approximately independent of the plasma parameters. 

Although the origin of the anomalous screening in the laboratory is unknown, if this 
phenomenon exists it should show up also in elastic collisions between nuclei, and thus it 
could affect the diffusion of nuclei through the plasma. This can be easily understood by 
observing that the anomalously large measured nuclear cross sections seem to indicate that 
nuclei are more transparent to each other. Actually, it is experimentally know that elastic 
cross sections of charghed nuclei are smaller than the Rutherford values at low energies 
(Huttel et al. 1985), however measurements of interst to stellar physics are so far not 
available. 

Again with a phenomenological approach, we know from nuclear fusion measurements 
that the Coulomb barrier is lowered by an amount Z7/ a &. Let us assume that a similar process 
occurs in the plasma, so that the effective interaction between two nuclei with charges Z\ 
and Z 2 at distance r is given by: 

V(r) = Z 1 Z 2 e 2 /r - U ph (14) 

The relative nuclear motion is thus the same as in a pure Coulomb potential, the collision 
energy being increased by a quantity U p i'. 

E -> E eff = E + U pl (15) 

The Coulomb cross section being cr Co „(E)= (ZiZ 2 e 2 / E) 2 , this means that in the plasma 
one has: 

crcouAE) = «Cou(E + U pl ) = (|x7T ) 2 - (16) 

For a simple estimate of this effect we note that the diffusion coefficient D are inversely 
proportional to the Coulomb cross section. By replacing E— > kT in eq. (16) one gets: 

D(U pl ) « D{0)[1 + . (17) 

A more precise estimate, properly taking into account the energy distribution for the colliding 
particle and still neglecting variation of the Coulomb logarithm, gives: 

/»oo e~ x x 2 

D(U pl ) = m I I <* % + W2 (18) 

The result of the numerical evaluation of the integral is well approximated by the following 
expression: 

D(U pl ) = D(0)[l + 3 I ^+ 1 -( 1 ^) 2 ] . (19) 
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This means that, in the region of interest to us [kT ss200 eV), an effective screening po- 
tential U p i of few hundred eV corresponds to diffusion coefficients increased by an order of 
magnitude. 

We do not attempt, here, to calculate all the diffusion coefficients relevant for the study 
of Li and Be. For this goal, at least one should have an estimate of the screening potential 
for the collisions of these ions with both H and He nuclei. It suffices here to remark that 
screening potentials of the same order as those measured in the laboratory can significantly 
enhance the rate of diffusion processes. 

4. LI AND BE IN HYADES LOW MAIN SEQUENCE STARS 

As already mentioned, surface abundances of light elements in low mass stars are the 
result of the concurrent action of two different physical processes: nuclear burning at the 
bottom of the convective envelope and atomic diffusion. It is generally difficult to discrimi- 
nate between the different contributions; however, when stellar ages are short enough with 
respect to the diffusion timescale, we have the opportunity to study the burning reaction 
rates for different values of the temperature at the bottom of the convective envelope in stars 
with different masses. This is the case of the Hyades, for which rather accurate measure- 
ments of Li and Be abundances are available (Boesgaard & Tripicco 1986; Duncan & Jones 
1983; Cayrel et al 1984; Thorburn et al. 1993; Boesgaard, Heacox and Conti 1977; Garcia 
Lopez et al. 1994). As is well known , although in the pre-main sequence the temperature at 
the bottom of the convective envelope reaches values larger then T& =3.5, the stellar lifetime 
is so short (less than 50 Myr) that classical stellar models cannot account for the observed 
Li depletion in the Hyades unevolved main sequence stars with T e ff < 6000 K. 

In Figs. 3-4 we report the surface abundance evolution of Li and Be for a 1M star with 
solar chemical composition calculated for various values of the effective screening potential. 

These results have been obtained by means of FRANEC (see Appendix A), neglecting 
diffusion. For all calculations presented in this section we use Z=0.02, Y=0.29 and the 
mixing length is a = 2.25. 

For Upf Ll =0 Li depletion is negligible, in agreement with previous well known results 
(Promt & Micheaud 1989, D'Antona & Mazzitelli 1984) but in contrast to observed abun- 
dances. Li abundance depletion in agreement with the observational result (i.e. Ay^ 1) 
can be obtained for Ufr u « (500 - 700) eV. 

Note that most of the variation of the Li abundance occurs in the pre-main sequence 
phase, a really negligible variation taking place between the end of pre-main sequence and the 
estimated age of this cluster {tjjy = 0.8 ± 0.2 Gy, see e.g. Castellani, ChiefR and Straniero 
1992), as the bottom of the convective region moves towards more external and cooler 
regions. In this context, Hyades abundances are representative of the pre-main sequences 
nucleosynthesis . 

An ad hoc effective screening potential could obviously reduce the amount of surface Li 
to the observed value in a star of a given mass. However it is not clear that the same value 
of Upf Ll can be adequate for stars of different masses. That this is the case is shown in Fig 
5, where we present the Li abundance for low mass Hyades as a function of their effective 
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temperature. A value of 1 ~ 700eV reproduces the trend of the observational data. It 
is worth noticing that for Upf L *=0 our results are again similar to those of other authors. 

Let us remark that the temperature at the bottom of the convective envelope depends, 
for a given stellar mass, on the assumed chemical composition. A larger metallicity could 
imply a higher temperature at the bottom (for a fixed value of the He abundance). We 
do not attempt here a detailed discussion of this point, and we only note that the value of 
Upi~ L * we just derived using Z=0.02 is, strictly speaking, an upper bound. 

No significant Be depletion is observed for Hyades stars with low effective temperatures 
(see Garcia Lopez et al. 1994). In this case, available data can be used to set an upper limit 
to the effective screening potential for p+Be reaction. Fig. 6 indicates that the observed Be 
abundances allow U^ Be < 1700 eV . 

5. FROM THE ZAMS TO THE PRESENT SUN 

Within a phenomenological approach, without resorting to detailed evolutionary calcula- 
tions, we study now the time evolution of the Li and Be abundances in the solar photosphere 
during the main sequence starting from Zero Age Main Sequence (ZAMS). Firstly we will 
write down the general equations for the evolution of heavy element abundances in the 
convective zone, showing that they are essentially determined from the properties of the in- 
nermost convective layer, essentially the temperature Tf,, the distance from the solar center 
Rb, the density p^, the pressure and temperature logarithmic derivatives, dlnP/ d(R/ R & ) 
and dlnT/d(R/R @ ). 

In the spirit expressed at the beginning of this section, we will assume that all these 
quantities change smoothly along the main sequence, their time dependence being of the 
form 

O a = O a@ [l + /3 a (t-t @ )/t @ ] (20) 

For the slopes j3 a we will use results of our evolutionary code including diffusion, see column 
b of Table 3, and the values &a at the bottom of the convective zone of the present Sun 
are determined from observational (helioseismological) data whenever possible, otherwise we 
will use results which are common to several solar model calculations, all this information 
being summarized in the Appendix B. 

In practice, to a very good approximation one could assume that all quantities, but the 
temperature T&, are constant. Since the nuclear burning rates vary strongly with tempera- 
ture, the time dependence of T& cannot be ignored and a few per cent decrease from ZAMS to 
the present Sun is important. We will take the temperature at the bottom of the convective 
layer in the present Sun, T&q, as a parameter in the range: 

T b@ = (2.1 - 2.3). (21) 

This interval is the result of inversion of helioseismological data (Christensen-Dalsgaard, 
Gough and Thompson 1991) and is consistent with evolutionary calculations where diffusion 
is considered, see again the Appendix B. 
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On these grounds we will calculate Li and Be abundances in a 1M & star of solar chemical 
composition along the main sequence. 



5.1 The basic equations 

Due to the short characteristic times of the convection processes, the abundances 
(yi=12-\-log w (Ni / Nh) and i=Li or Be) of each element in the photosphere are the same as 
in any point of the convective zone. The time evolution of can be represented as: 

^ = (log w e)^^ = (log w e)(- < X p p f > +v l p l S/M l ) i = 7 Li, 9 Be (22) 

where Mi is the mass of the i-th element contained in the convective zone, and the two terms 
on the right hand side, which we are going to define and to discuss in the following, represent 
respectively the contributions to the i-th element depletion due to the nuclear burning in 
the convective zone and due to the outflow towards the inner radiative region. 

Coming to the first term, < \ P J^ 1 > is the nuclear burning rate averaged over the convec- 
tive zone 



<Xpl >~ JdM ' (23) 



where the integral is over the convective zone and the nuclear burning rates \ P J^ 1 are given 
by equations (4) and (8). 

Due to the strong temperature dependence, nuclear burning occurs most likely at the 
bottom of the convective zone, and its average rate is essentially related to the properties 
of the innermost layer. By numerical experiments (performed at different solar ages and by 
varying the position of the bottom of the convective zone) we find, to an accuracy of about 
10%, 



i„ ,u: 



p+i 
pi 



bottom 8 H b y kT b J y ' 



and in the following we will fix = Xb = 0.7, Pb&= 0.2 g cm -3 . 

Coming to the last term in equation (22), Vi and pi are respectively the diffusion velocity 
and the density of the i-th element calculated at the bottom of the convective layer and 
S=47ri?^ is the surface defining the border between the convective and the radiative region. 
We will take: 

R b& = 0.71R & (25) 
Clearly the factor pi/Mi does not depend on the element one is considering, and one has: 

Pi I Mi = Pb /M conv (26) 
where M conv is the mass contained in the convective zone; we will take M conv& = O.O25M . 
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Diffusion of the i-th element is driven by pressure gradients, thermal gradients and con- 
centration gradients. At the bottom of the convective region the concentration gradients of 
H and He are negligible, and we only take into account Li (Be) concentration gradients, as 
these elements can be burnt in radiative regions which are closeby. We write thus for the 
diffusion velocity Vi = Vi • R: 

„ r . dlnCi . dlnP 4 dlnT , . „. 

"■' = + " + iT m)' (27) 

where the scale factor D Q can be written as: 

D = 1.2 • 10" 11 T 5/2 /p b cm/sec (28) 

the density being measured in c.g.s units 

The coefficients Ap,Ax and Ac are weakly dependent on temperature and density, at 
least in the region of interest to us, and we will take them as constant. The values for 
parameters representative of the present solar convective basis have been calculated following 
Thoul et al. (1994) and are shown in Table 4. 

One notes that the Li and Be coefficients are similar and that for both elements the 
concentration coefficient is typically a factor 10 smaller than that of pressure. It has to be 
remarked, for future applications, that this is not an accident, but rather it is a consequence 
of the general relation between mobility and diffusion coefficients (see for example Landau 
& Lifshitz, §12 of Physical Kinetics and §59 of Fluid Mechanics). Pressure and thermal 
coefficients are comparable; however, as is well known, in the Sun thermal diffusion is less 
important than the pressure term, the thermal gradient being significantly smaller than that 
of pressure. 

As discussed in section 3, the actual diffusion coefficients could differ from those calcu- 
lated above due to the effects of screening, or to some turbulent process. For this reason, 
we will introduce a multiplicative parameter n in front of equation (28): 

D ^D = V D , (29) 

which we assume to be the same for Li and Be and which we will fix later from the observed 
Be depletion. Coming to the pressure and thermal gradients in equation (27), we note (see 
Appendix B) that for the present Sun all calculations give similar values at the bottom of 
the convective zone, close to: 



dlnP 



d(R/R @ ) 



-12 (30a) 



dlnT 
d(R/R & ) 



(30b) 



We can approximately estimate the concentration gradient, by assuming that at any 
time the concentration below the bottom of the convective zone is determined by the nuclear 
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burning rate only, i.e. as a first approximation we neglect diffusion since this process is slow 
in comparison with the others. In this way one gets: 

C i (t) = C i (0)ex P [-£dt'\ p + i ] (31) 

where time corresponds to ZAMS. It follows that dlnCi/d(R/ R & ) can be expressed in 
terms of the burning rate and temperature gradient at the bottom of the convective zone: 

In the approximation that the properties of the basis of the convective zone are unchanged 
during the main sequence, the above equation simplifies to: 

dlnC ' -tKr^-Krimj^^. (33) 



d(R/R @ ) pi i 1 P i , J(rfi2/i2 ) 



5.2 Results 

By integrating numerically equation (22) from ZAMS to the present Sun, we find the 
results shown in Figs. 7-10, which deserve several comments. 

i) As expected, see Fig. 7 dashed curve, no significant Li depletion during the main sequence 
occurs, for the conventional (U^ Ll = 0) nuclear burning rates and for standard (7/ = 1) 
diffusion coefficients. 

ii) For the same value which reproduces the observed Li abundance in 1M Hyades stars, 
t/£ +Li 700 eV, Li is burnt durning the main sequence but not enough to account for the 
present solar value. Standard diffusion terms (7/ = 1) contribute little to Li depletion, see 
Fig. 7. 

iii) For large enough diffusion coefficients (7/ = 5), the Li depletion increases, in particular 
we note that the contribution of the concentration driven diffusion is relevant, see Fig. 8. 

iv) For standard diffusion coefficients, the Be abundance is unchanged from the Hyades age 
to the present Sun, even for very large values of U^ Be , see Fig. 9 full curves. In other 
words, Be cannot be burnt. On the other hand, for 7/ 5 (dashed curves in Fig. 9) diffusion 
is effective to bring below the convective zone an amount in agreement with observational 
data. We remark that 7/ 5 corresponds to screening potentials of the order of few hundred 
eV, see section 3. We thus empirically choose 7/ = 5, for both Be and Li. 

v) It is worth observing, see Fig. 10, that for 77 ~ 5 and 

Ki Ll ~ 700 eV the Li de P letion 

during the main sequence is Ay^ 1 in agreement with observational data. 

We have also evaluated the Li surface abundance evolution in a 1M star directly with 
a recent version of the FRANEC code where diffusion of He and heavy elements (C,N,0, Li 
and Be) is included. Diffusion coefficients have been calculated as in Thoule et al. (1994), a 
multiplicative factor 7/ = 5 being included for Li and Be. The results are presented in Fig. 



13 



11, where one sees that again for Ujf « 600 - 700 eV we are able to reproduce all the 
observational data from meteorites to the present Sun. 

We remark that this range of V^ L% corresponds to that needed for explaining the Li 
abundance in low main sequence Hyades, see section. 4 and it is not far from the value 
measured in the laboratory Uf+ Li « 400 eV. 

6. CONCLUSIONS AND OUTLOOK 

Briefly, the main result of this paper is the following: if one assumes that nuclei in the 
stellar plasma are more transparent to each other, as suggested by laboratory experiments 
in the 10-100 KeV region, one can understand light (Li and Be) elements abundances both 
in the Sun and in lower main sequence Hyades. 

Quantitative results are summarized in Fig. 5, where we compare the calculated and 
the observed Li abundances in low main sequence Hyades, and in Fig. 11, where we show 
the time evolution of Li abundance in 1M stars, calculated by using our FRANEC code 
(element diffusion being included), together with the observational values. 

In short, for effective plasma screening potentials of the order of few hundred eV the 
Li burning temperature is reduced so that it can be burnt at the bottom of the convective 
layer and at the same time Be diffusion is enhanced so that it can be hidden below the solar 
convective zone. 

We admit that the proposed connection between the laboratory and stellar interiors is 
highly hypothetical and that so far there is no explanation for the enhancement of astro- 
physical S-factors at low energies, nevertheless this enhancement is now the result of several 
experiments and the values we derive for the effective plasma screening potential are rela- 
tively similar to those measured in the laboratory. 

We believe that a series of new experiments and observations is necessary for test- 
ing/disproving our approach: 

i) the cross sections of p+ 7 Li — > ck + ck, p+ 9 Be — > a + 6 Li , p+ 9 Be — > d-\-2a reactions should 
be measured at still lower energies, so that a precise determination of Uf^ Ll and Uf^ Be can 
be obtained; 

ii) elastic (transport) cross sections for p+ 7 Li and p+ 9 Be collisions should be measured in 
an energy range of astrophysical interest: does the enhanced transparency suggested by 
inelastic reactions also show up in the elastic channel? 

iii) precise measurements of Li and Be abundances in intermediate age and old open clusters 
would obviously be interesting, to compare the time dependence we propose with richer data 
set. 
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APPENDIX A: THE EVOLUTIONARY CODE 



The stellar models presented in this paper have been computed by means of the latest 
version of the FRANEC (Frascati Raphson Newton Evolutionary Code), described by Chieffi 
& Straniero (1989; CS89). Let us recall the relevant physical inputs and computational 
procedures and comment on some recent improvement. 

This version of the code differs from that described in CS89, mainly because the equations 
of chemical evolution and stellar structure are simultaneously solved. In addition a full 
network that includes all stable isotopes up to 64 Ni is explicitly included. 

The equation of state takes into account the quantum-relativistic effects for the electron 
component of the stellar plasma and the electrostatic interaction (Straniero 1988). 

Tables of radiative opacity coefficients are derived form Iglesias, Roger and Wilson (1992) 
for T> 10 4 K and from Kurucz (1991) at lower temperature. The heavy element solar 
mixture of Grevesse (1991) has been adopted (except for low-temperature opacity for which 
a slightly different mixture has been used by Kurucz, namely Anders & Grevesse 1989). 

Nuclear reaction rates are generally taken from Caughlan & Fowler (1988), whereas 
for the 7 Li(p,ct) reaction we have adopted the more recent cross section measurements by 
Engstler et al. 1992. Electron screening is derived form Graboske et al. 1973 (see also 
DeWitt, Graboske and Cooper 1973). 

As usual, the mixing length has been calibrated by comparing the theoretical radius of 
the standard solar models (SSM) with the observed one (see Chieffi, Straniero and Salaris, 
1995). In such a way we have obtained a = 1/ H p =2.25 for the no-diffusion SSM and 
a = 2.40 when diffusion of H, He, C, N, 0, Li and Be are taken into account following 
Thoule et al. 1994. The other relevant features of our SSM are reported in Table 5. 

The allowed variations of r, L, P, T, M(r) and temperature gradient between two adjacent 
mesh points should not be in excess of some prefixed values, namely: Sr/r=0.1, SL/ L =0.01, 
6P/P=0M, £T/T=0.02, SM/M=0.01 and £VT/VT=0.1. As a consequence, the typical 
number of mesh points required for an integration for a PMS or a MS models ranges between 
500 and 600. About 600 time-steps are required for a SSM, 400 of which for the PMS. 

APPENDIX B: THE BOTTOM OF THE SOLAR CONVECTIVE ZONE 

We present in Table 5 a summary of information about the bottom of the convective zone 
in the present Sun, as obtained from helioseismological observations and from the results 
of several standard solar model calculations. Concerning our results (F), these have been 
obtained by using the FRANEC code with and without diffusion (see Appendix A). 

From inversion of the solar data (Christensen-Dalsgaard et al. 1991) it has been possible 
to derive the depth of the convective zone (<4= 1- Rb/R®) and the sound speed, q,. The 
temperature is then estimated from the sound speed, by using the perfect gas law and 
assuming the He mass fraction to be in the range Y=0.23-0.29: 

T b = (2.2 ± 0.1) 10 6 K (34) 

Concerning standard solar model calculations, when T& (q,) was not given by the authors, 
it was calculated by us from the value of q (T&) again by using the fully ionized perfect gas 
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law. 

As is well known (Bahcall & Pinsonneault 1992; Christensen-Dalsgaard, Proffit and 
Thompson 1993; Promt 1994), in models without diffusion the convective zone is thinner 
and the basis is cooler than given from helioseismological data. On the other hand, models 
where diffusion is taken into account look closer to the real Sun. 

Coming to the quantities which are relevant for diffusion, we note that the gradients 
dlnP I d(R/ Rq) and dlnT / d(R / R & ) are similar (to the level of 10% or better) among the 
different models. 

We define a P-T diffusion rate parameter as: 

= log w e v lPT piS/Mi , (35) 

PT 

where Vipj is the diffusion velocity due to pressure and thermal gradient, see equation (27). 
It is important that (^)px does not change by more than 20% when different solar models 
are used. 

The time dependence (along the main sequence) of the physical parameters, O a at the 
bottom of the convective zone has been investigated numerically within the FRANEC code 
(with and without diffusion). As expected all quantities vary smoothly and weakly and are 
well approximated by a linear function of time: 

O a = O a@ [l + /3 a (t-t @ )/t @ ] . (36) 

In Table 3 the values of j3 a are shown. We remark that nuclear burning rates are very 
sensitive to T&. Thus even a 10% variation along the main sequence can be significant. 

On the other hand we note that the diffusion rate parameters (^f)pr are constant in 
time within 10% or better. From Table 5 one thus derives that the photospheric Li and Be 
depletions during the solar main sequence are, for rj = 5: 

AyfJ (MS) « Ay^(MS) « 0.06 . (37) 
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TABLES 



TABLE I. Data on Be abundances (y Be = log w (N Be /N H ) + 12). 

VBe Ref- 

Meteorites 1.42± 0.04 Anders & Grevesse (1989) 

1M star in the Hyades 0.9±0.3 Garcia Lopez et al. (1994) 

Solar photosphere 1.15± 0.10 Anders & Grevesse (1989) 



TABLE 2. The experimentally determined effective screening potential, Ui ab , together with the 
theoretical value in the adiabatic limit, U e i, calculated following equation (13). 



reaction 


U lab [eV] 


U d [eV] 


Ref. 




atomic target molecular target 








186 ± 9 
470±150 
300±280 
380±250 



123±9 
440±150 
300±160 
330±120 

0-2000*? 

430±80 



^110 
;180 
;180 
;180 

;260 

;345 



(a) 
(b) 
(b) 
(b) 

(c) 

(d) 



(a) Prati et al. 1994, (b)=Engstler et al 1992, (c)=Tombrello & Sierk 1989, (d)=Angulo et al 1993. 
* value estimated by us. 



TABLE 3. Time dependence for the physical parameters, O a , of the Sun during the main 
sequence phase, at the bottom of the convective zone, We present the values of the slopes {3 a in 
the linear approximation O a = Oq,q[1 + f3 a (t — *©)/*©] as calculated through FRANEC without 
(a) and with (b) diffusion. 

(a) ~W 
T b -0.10 -0.07 

p b -0.75 -0.70 

R b 0.1 0.1 

M conv -0.25 -0.13 

dlnP/d(R/R @ ) -0.15 -0.16 

dlnT/d(R/R & ) -0.12 -0.18 

{dy Ll /dt) PT 0.06 0.01 
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TABLE 4. Diffusion coefficients at the bottom of the solar convective zone, calculated for 
T b = 2.2 • 10 6 K, p = 0.2 g/cm 3 and X b = 0.7, following Thoule et al 1994. 

element A p At A c 

^Li 0^9 IT2 -0.1 

9 Be 0.8 1.3 -0.06 



TABLE 5. Comparison among different standard solar models without and with diffusion; the 
helioseismological results (H) are also shown. The labels correspond to: H=Cristensen-Dalsgaard, 
Gough and Thompson (1991); CPT=Christensen-Dalsgaard, Proffitt and Thompson (1993); P = 
Proffitt (1994); TCL=Turck-Chieze & Lopes (1993); BP=Bahcall & Pinsonneault (1992). F is for 
the FRANEC code we are using. 





H 


F 


CPT 


TCL 


P 


BP 


CPT 


F 


P 


P 


BP 


Fixed parameter 




Z 


Z 


Z/X 


Z/X 


Z/X 


Z 


Z 


Z/X 


Z/X 


Z/X 


Diffusion 




no 


no 


no 


no 


no 


He 


He,Z 


He 


He,Z 


He 


L [lO 33 erg/ S ] 




3.83 


3.85 


3.85 


3.85 


3.86 


3.85 


3.83 


3.85 


3.85 


3.86 


R o [10 10 cm] 




6.96 


6.96 


6.96 


6.96 


6.96 


6.96 


6.96 


6.96 


6.96 


6.96 


t@[Gy] 




4.6 


4.6 


4.5 


4.63 


4.6 


4.6 


4.6 


4.63 


4.63 


4.6 






2.910 


2.798 


2.714 


2.729 


2.716 


2.777 


2.890 


2.740 


2.803 


2.727 


Z iniz x 1Q 2 




2. 


2. 


1.770 


1.907 


1.895 


2. 


2. 


1.988 


2.127 


1.958 


(Z/X) im ' z X 10 2 




2.903 


2.856 


2.490 


2.694 


2.671 


2.848 


2.894 


2.802 


3.045 


2.766 


ysurf x 1Q 




2.910 


2.798 


2.714 


2.724 


2.716 


2.494 


2.552 


2.456 


2.514 


2.466 


Z surf x 1Q 2 




2. 


2. 


1.770 


1.907 


1.895 


2. 


1.849 


1.979 


1.964 


1.958 


(Z/X) 3 "^ X 10 2 




2.903 


2.856 


2.490 


2.694 


2.671 


2.737 


2.546 


2.694 


2.694 


2.668 


(R b /R Q ) x 10 


7.13±0.03 


7.29 


7.22 


7.25 


7.25 


7.21 


7.08 


7.13 


7.103 


7.115 


7.07 


Cf,[l0 7 cm/s] 


2.23±0.02 


2.14 


2.23 


2.23 




2.19* 


2.23 


2.23 






2.25* 


T b [W 6 K] 


2.21±0.10 


2.09 


2.22* 


2.11 




2.13 


2.22* 


2.21 






2.26 


p b [g/cm 3 ] 




0.154 








0.167 




0.179 






0.197 


(M co „„/M ) X 10 2 




2.03 






2.06 


2.16 




2.35 


2.48 


2.48 


2.54 


dlnP/d(R/RQ) 




-13 








-12 




-13 






-12 


dlnT/d(R/R Q ) 




- 4.9 








- 4.6 




- 4.7 






- 5.1 


(dvLi/d^pTlGy- 1 ] 




- 0.014t 








- 0.013t 




- 0.013t 






- 0.013t 


{dy Be ldt) PT [Gy-^] 




- 0.014t 








- 0.012t 




- 0.012t 






- 0.012t 



* caicuiated by us assuming fuiiy ionized perfect gas law. 

f calculated by us with r) = 1 and diffusion coefficients of Table 4. 
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FIGURES 



FIG. 1. Observed Li abundances, yn = log^Nn/Nn) + 12, in 1M stars as a function of 
age. For the open clusters, the error bars are indicative of the spread in the measured values 
and/or of the uncertainties on the measured T e ff- Also shown are the measured value in the 
solar photosphere and the meteoritic value. The straight line is a linear best fit to all data, but 
the meteoritic. On the right side we also show the estimated pre-main sequence (PMS) and main 
seqence (MS) Li depletion (Ay^i) for the Sun. 

FIG. 2. The reactivity Na < av > for collisions of Li and Be with protons as a function of 
temperature. For 9 Be we use the expression given by Caughlan & Fowler 1988; for 6 Li and 7 Li we 
take the more recent results from Engstler et al. 1992. 

FIG. 3. Time dependence of the Li abundance, yn = log^Nn/Nn) + 12, in a 1M stars 
with solar chemical composition, calculated for different values of the effective plasma screening 
potential U^ L% . Diffusion effects are not taken into account in this calculation. Also shown are 
the meteoritic value, taken as the initial stellar abundance, and the observational result for 1M 
stars in the Hyades. 

FIG. 4. Time dependence of the Be abundance, ys e = i°9(^Be/^H) + 12, in a 1M stars 
with solar chemical composition, calculated for different values of the effective plasma screening 
potential U^ Be . Diffusion effects are not taken into account in this calculation. Also shown is the 
meteoritic value, taken as the initial stellar abundance. 

FIG. 5. Li abundance, yn, for the Hyades dwarfs as a function of the effective temperature, 
T e ff- We show data from Boesgaard & Tripicco (1986), circles; from Duncan & Jones (1983), 
diamonds; from Cayrel et al. (1984), squares, and from Thorburn et al. (1993), diagonal crosses. 
The cross on the top right is indicative of the observational errors. The three curves represent the 
calculated behaviour for the values of U^ Ll [eV] shown in the figure. All curves are normalized 
to yield (approximately) the observed Li abundance at T e j^=6300. 

FIG. 6. Be abundance, yse, f° r the Hyades dwarfs as a function of the effective temperature, 
T e ff. We show data from Garcia Lopez et al. (1994), filled circles, and from Boesgaard, Heacox 
and Conti (1977), open circles. The three curves represent the calculated behaviour for the values 
of Upf Be [eV] shown in the figure. All curves are normalized to yield (approximately) the observed 
Be abundance at T e j^=6300. 

FIG. 7. Time dependence of the Li abundance from ZAMS to the present Sun, calculated 
assuming T{, = 2.2 • 10 6 K, U^ Ll = 700 eV and 77=1. We show the contributions due to nuclear 
burning (BUR), pressure and thermal diffusion (PT), concentration driven diffusion (CON), and 
the total (TOT). We also show the result of the standard (U^ Ll = and 77 = 1) calculation 
(dashed line, STA). 
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FIG. 8. Time dependence of the Li abundance from ZAMS to the present Sun, calculated 
assuming T b& = 2.2 • 10 6 K, C/^ +L ' = 700 eV and 77=5. Same notation as in Fig. 7 

FIG. 9. Variation of the Be abundance between the Sun age and the ZAMS as a function of 
Upf Be , for different values of T&q, calculated: 

a) for 77= 1 (full lines), 

b) for 77= 5 (dashed lines) 

FIG. 10. Variation of the Li abundance between the Sun age and the ZAMS a function of 
Upf Be , for different values of T&q, calculated: 

a) for 77= 1 (full lines), 

b) for 77= 5 (dashed lines) 

FIG. 11. Time dependence of the Li abundance, yn = log^Nn/Nn) + 12, in a 1M star with 
solar chemical composition, calculated for different values of the effective plasma screening potential 
Upi~ L * m e V- Diffusion effects are taken into account, see text. The multiplicative factor for Li and 
Be diffusion is set v)Li,Be = ^- Also shown are the observational data. 
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